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We study the low-temperature properties of a classical Heisenberg model with site-random in- 
terlayer couplings on the cubic lattice. This model is introduced as a simplified effective model 
of Sr(Fei_2,Mn2,)02, which was recently synthesized. In this material, when x — 0.3, (mrn) and 
(ttttO) mixed ordering is observed by neutron diffraction measurements. By Monte Carlo simula- 
tions, we find an exotic bulk spin structure that explains the experimentally obtained results. We 
name this spin structure the "random fan-out state". The mean-field calculations provide an in- 
tuitive understanding of this phase being induced by the site-random interlayer couplings. Since 
Rietveld analysis assuming the random fan-out state agrees well with the neutron diffraction pattern 
of Sr(Feo.7Mno.3)02, we conclude that the random fan-out state is reasonable for the spin-ordering 
pattern of Sr(Feo.7Mno.3)02 at the low-temperature phase. 
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I. INTRODUCTION 

Disorder in materials induces the frustration effect and 
makes it difficult for cooperative phenomena to occur [ll- 
0]. For example, in random magnets, conventional mag- 
netic orders such as ferromagnetic (FM) and Neel orders 
are suppressed by the frustration that arises from disor- 
der. Consequently, the ordering temperature of random 
systems is lower than that of pure systems and the or- 
dered phase disappears due to random interactions be- 
tween magnetic moments in some cases. Systems with 
disorder also have many features that are not observed 
in pure systems. In some random magnetic systems, 
interesting phases and transitions have been observed, 
for example, s pin gla ss (7|-[l0j. mixed ordering , and 

oblique phase [17'->19]. To find phases having exotic static 
and dynamic behavior, synthesis and characterization of 
random magnets have been attracted attention in mate- 
rials science for a long time. 

From the viewpoint of statistical physics, it is impor- 
tant to explore how disorder in magnetic systems af- 
fects phase transitions and magnetic orders. A theo- 
retical model of random magnets constructed by ran- 
dom substitution of magnetic ions is a site-random 
modelfH, [H, [l^ l20l - [26| . Regarding a binary mag- 
netic alloy, we can describe a site-random model with two 
types of magnetic ions, A and B. In this model, A and 
B ions are placed randomly on the lattice. The interac- 
tions of each bond depend on the combination of ions. 
For example, the interaction between A ions is antiferro- 
magnetic (AF), whereas other interactions are FM. The 
Heisenberg model using this decision rule of interactions 
for all directions on a three-dimensional lattice has been 
investigated by Monte Carlo simulations[12, 13, 15, 16]. 
This model can be considered as a fundamental model of 
isotropic materials. However, we cannot adapt this model 



in the present case because origins of interaction depend 
on the direction in the infinite-layer structure treated in 
this paper. In this case, the emergence of disorder of 
interaction should depend on the direction. In isotropic 
random systems, the key element is the concentration 
density of each ion. In anisotropic random systems, in 
addition to that, the spatial distribution of disorder is 
also important. Then it is an interesting issue to inves- 
tigate phase transition and spin structure in anisotropic 
random systems. Since recent developments in synthetic 
methodology enable us to design the spatial structure of 
materials, the above-mentioned issue has become increas- 
ingly important in materials science. 

Recently, the infinite-layer iron oxide SrFe02 was syn- 
thesized by hydride reduction of SrFeOs, which has 
the perovskite structure f27j|. From the crystal struc- 
ture of SrFe02, the nearest-neighbor spin interaction 
J in a plane is due to a Fe-O-Fe superexchange in- 
teraction, and the interlayer interaction J' is due to 
direct through-space overlap between Fe ions. First- 
principle calculation s l28l [29j and inelastic neutron scat- 
tering measurements [30j consistently revealed that both 
interactions are AF with J ~ 3 meV and J' ~ 1 meV. 
Owing to these strong AF interactions, SrFe02 exhibits 
(tttttt) order (see Fig. [1] (a)) with a relatively high Neel 
temperature Tn — 473 K. Seinberg et al. studied the 
effect of substituting Mn^+ ions (d^) for Fe^+ ions (d^) 
at room temperature (RT)|31J. RT neutron diffraction 
and Mossbauer spectroscopy studies on Sr(Fei_a;Mna;)02 
{x = 0.0, 0.1, 0.2, 0.3) showed that Mn substitution at 
Fe sites substantially destabilizes (tttttt) order with in- 
creasing x, and Sr(Feo.7Mno.3)02 becomes paraniagnetic 
at RT[3iJ. From the Goodenough-Kanamori rule ^,[3^, 
the Fe-O-Mn superexchange intralayer interaction must 
be AF. Thus, the destabilization of (tttttt) order suggests 
that the Fe-Mn interlayer interaction is FM and com- 
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FIG. 1: (Color online) (a) Schematic spin structure of (mrn) 
order, (b) Schematic spin structure of (ttttO) order. 



petes with the AF Fe-Fe interlayer interaction J'[3l|. 
From the above facts, the material Sr(Fei_a;Mna;)02 is 
an anisotropic random system. It is possible that an un- 
usual phase transition and a spin structure that comes 
from interlayer random couplings emerge at low temper- 
ature. 

This paper presents low-temperature neutron diffrac- 
tion results on Sr(Feo.7Mno.3)02, as well as a theoret- 
ical interpretation of those results. We find that two 
distinct wave vectors, namely, (tttttt) and (ttttO), de- 
velop simultaneously in Sr(Feo.7Mno.3)02 at low temper- 
ature. The aims of our study are to clarify the spin- 
ordering pattern and to explain the mechanism whereby 
the (tttttt) and (ttttO) wave vectors coexist. In Sec. II, we 
present the experimental results on the random magnet 
Sr(Fei_a;Mna;)02. In Sec. Ill, we introduce a classical 
Heisenberg model with site-random interlayer couplings 
as a simplified effective model of Sr(Fei_j;Mn2;)02. In 
Sec. IV, we show the Monte Carlo simulation results 
of finite-temperature properties of the model introduced 
in Sec. III. Drawing the phase diagram of temperature 
versus Mn ion concentration, we find the mixed phase 
which is characterized by wave vectors (tttttt) and (ttttO). 
We discuss the spin-ordering pattern of this phase and 
the universality class of each phase transition. In Sec. 
V, to consider the emergence mechanism of the mixed 
phase observed in Monte Carlo simulations, we investi- 
gate the effect of random interlayer couplings by mean- 
field calculations. In Sec. VI, we discuss whether the 
spin-ordering pattern is robust with respect to the de- 
tails of the model, such as decision rules of interactions 
and easy-plane anisotropy. We also show the results of 
Rietveld analysis for Sr(Feo.7Mno.3)02. Section VII sum- 
marizes the paper and our main conclusions. 



II. EXPERIMENTAL RESULTS 

The random magnet Sr(Fei_3;Mna;)02 (x = O.I, 0.2, 
0.3) was recently synthesized [3l|. The crystal structure 
of Sr(Fei_a;Mna;)02 is shown in Fig. [5] As we pointed 



FIG. 2: (Color online) (a) Crystal structure of 
Sr(Fei_a;Mna;)02. Orange (large), blue (medium), and 
red (small) spheres indicate Mn, Fe, and O ions, respectively. 
This arrangement of Fe and Mn ions is only an illustrative 
example. In Sr(Fei_2,Mn2,)02, Fe and Mn ions are randomly 
placed on the tetragonal lattice. 



out, the a; = G.f and 0.2 samples exhibit (tttttt) order 
above RT. The transition temperatures of the two sub- 
stituted compounds have not been investigated, but the 
experimental results in Ref. |3l| indicate that Tn de- 
creases with increasing x. Powder neutron diffraction ex- 
periments at RT show a systematic decrease in magnetic 
reflection intensity; specifically, the relative intensity of 
the (1/2,1/2,1/2) reflection for x = 0.0, 0.1, and 0.2 is 
1:0.8:0.4. ^^Fe Mossbauer experiments at RT also reveal 
a systematic decrease of hyperfine fields: 402 kOe for x 
= 0.0, 330 kOe for x = 0.1, and 26.1 kOe for x = 0.2. 
When X — 0.3, the system is in a disordered state at RT. 
Thus, we performed neutron diffraction measurements at 
low temperature. 

Powder neutron diffraction experiments on 
Sr(Feo.7Mno.3)02 were carried out on a Kinken powder 
diffractometer with multicounters for HERMES (High 
Efficiency and high Resolution MEasurementS) at the 
Institute for Materials Research, Tohoku University, in- 
stalled at a guide hall of the JRR-3 reactor of the Japan 
Atomic Energy Agency (JAEA), Tokai 34]. The incident 
neutrons were monochromatized to 1.8204 A by the 331 
reflection of bent Ge crystal. A 12'-blank-sample-18' 
coUimation was employed. A polycrystalline sample of 3 
g was placed into an He-filled vanadium cylinder. The 
sample temperature was controlled from 10 K to 273 
K. The data were collected with a step-scan procedure 
using 150 neutron detectors in a 29 range from 3° to 
153° with a step width of 0.1°. 

Figure [3] shows the neutron diffraction pattern for 
Sr(Feo.7Mno,3)02 at 10 K. As expected from the previ- 
ous report [3l|, the diffraction profile has magnetic peaks 
such as (1/2,1/2,1/2) and (1/2,3/2,1/2) characterized by 
a wave vector (tttttt), together with a nuclear reflection 
with a tetragonal PA/mmm cell (a = 4.001 A, c = 3.441 
A). However, we unexpectedly observed magnetic peaks 
such as (1/2,1/2,0) and (1/2,1/2,1) characterized by a 
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FIG. 3: (Color online) Powder neutron diffraction pattern 
of Sr(Feo.7Mno.3)02 at 10 K. Purple and blue indices cor- 
respond to magnetic reflections of (ttttO) and (tttttt) orders, 
respectively. Triangles correspond to peaks from perovskite 
Sr(Feo.7Mno. 3)03-5, and asterisk corresponds to a peak from 
an unknown impurity. 
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FIG. 4: (Color online) Temperature dependences of 
(1/2,1/2,1/2) and (1/2,1/2,0) peaks corresponding, respec- 
tively, to (tttttt) and (ttttO) wave vectors for Sr(Feo.7Mno.3)02. 
Dotted curves are visual guides. The intensities are normal- 
ized by acquisition time. 



wave vector (tfttO) at the same time. Plots versus tem- 
perature of the (1/2,1/2,1/2) and (1/2,1/2,0) peaks cor- 
responding, respectively, to the (tttttt) and (ttttO) wave 
vectors are shown in Fig. U) Here, it is seen that the 
two peaks lose their intensities gradual with increasing 
temperature, and at and above 240 K, both peaks dis- 
appear and the intensities decrease to the background 
level. The pure SrFe02 system exhibits (tttttt) order 
only, and thus the appearance of the (ttttO) peak is 
caused by the effect of substitution on Mn ions. We can- 
not predict the sign of direct exchange interaction along 
the interlayer direction, because the direct exchange in- 
teraction strongly depends on the environment around 
it. In contrast to the interlayer interaction, from the 
Goodenough-Kanamorirule,35], Fe-O-Mn and Mn-O-Mn 
superexchange intralayer interactions must be antiferro- 
magnetic. Thus, the simultaneous appearance of (tttttt) 
and (ttttO) wave vectors indicates that interlayer FM in- 
teractions appear at Mn sites. 

A simple explanation of the observed simultaneous 
(tttttt) and (ttttO) wave vectors is phase separation into 
the (7TTTTT)-type spin structure (Fig.[T](a)) and the (ttttO)- 
type spin structure (Fig.[T](b)). The other scenario is the 
emergence of a bulk spin structure in which (tttttt) and 
(ttttO) wave vectors coexist. Since Sr(Fei_2;Mna;)02 is an 
anisotropic random system that has never been investi- 
gated theoretically, it is possible this type of bulk spin 
structure emerged in Sr(Fei_a;Mna;)02. Indeed, as we 
will see in the following theoretical argument, we find a 
novel type of magnetic structure that can explain the si- 
multaneous appearance of (tttttt) and (ttttO) wave vectors. 
From these considerations, in Sec. VI B, we discuss the 
spin-ordering pattern in Sr(Feo.7Mno.3)02 at low tem- 
perature. 



III. THEORETICAL MODEL 

To investigate whether there is a bulk spin struc- 
ture that explains the simultaneous appearance of 
(tttttt) and (ttttO) wave vectors without any phase sep- 
aration, we construct a simplified effective model of 
Sr(Fei_a;Mna;)02. We introduce a classical Heisenberg 
model with site-random interlayer couplings on the cu- 
bic lattice. This model includes the characteristic fea- 
tures of Sr(Fei_a;Mna;)02: (i) the system consists of two 
types of magnetic ions randomly placed on the lattice and 
(ii) competition between FM and AF interactions exists 
along only the interlayer direction. The model Hamilto- 
nian is given by 

H^J ^ s^-Sj + J' ^ aijS^-Sj, (1) 

intralayer interlayer 

:=-i{l + (e, + ej)-e»£j}, (2) 

where S; is the three-component vector spin of unit length 
at the i-th site on the cubic lattice, and Si is equal 
to —1 or -1-1 for A or B ions at the z-th site, respec- 
tively Here, A and B ions correspond to Fe and Mn 
ions in Sr(Fei_2:Mn2:)02. The first term of the Hamil- 
tonian denotes the uniform AF nearest-neighbor interac- 
tion (J > 0) in the afe-plane. The sign of the nearest- 
neighbor interaction in the second term depends on the 
arrangement of A and B ions along the interlayer di- 
rection (c-axis): the interactions between A ions are as- 
sumed to be AF where cjij = +1, and other interactions 
are assumed to be FM where cr.y = — 1. We assume that 
absolute values of the interactions are the same {J' — J) 
for simplicity. The interactions of this model are sum- 
marized in Table HI Under this rule of interactions, we 
assume that the interlayer interaction at B ion sites is 
always FM. However, there are other possible decision 
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rules of interactions where FM interactions exist at B 
ion sites. In Sec. VI, we discuss the low-temperature 
properties depending on the decision rule of interactions. 

We calculate average densities of FM and AF interac- 
tions along the interlayer direction depending on a B ion 
concentration, x. The probabilities of A-A, A-B, and B- 
B for nearest-neighbor ions along the interlayer direction 
are (1 — xY, 2x{\ — x), and x^, respectively. Thus, prob- 
abilities of FM and AF interactions are x^ + 2x{l — x) 
and (1 — x)^, respectively. The spatial average of signs 
of interlayer interactions, ct, is given by 

a = 2x^ - Ax + 1, (3) 

where a is equal to Q &i x ^ x* (2 - %/2)/2 ~ 0.2929. 
In the case of x < x* , the AF interaction is dominant 
(ct > 0) along the interlayer direction. In contrast, the 
FM interaction is dominant (ct < 0) for x > x* . 

In the site-random model, the signs of interactions 
along the interlayer direction are correlated, and thus 
the FM correlation between next-nearest layers (NNLs) 
emerges. In fact, the effective FM interaction between 
NNLs exists regardless of x in our model. The sign of 
the effective interaction between next-nearest-neighbor 
sites is defined by 

c^jT' := -CTijCTjfc, (4) 

where i and k denote next-nearest-neighbor sites along 
the interlayer direction through the site j. To show 
that effective interactions are always FM, we investigate 
the arrangements of three ions including two interactions 
along the interlayer direction. If the signs of both inter- 
actions are the same, that is, Uij — ajk, the effective 
interaction is FM (cr"^^' = — 1). In contrast, if one of 
the interactions is FM coupling and the other is AF cou- 
pling, that is, (Ty = ~'^j7c, the effective interaction is 
AF (cr™' = -1-1). For example, the effective interactions 
of arrangements A-A-A and A- A-B are FM and AF, re- 
spectively. We summarize the effective interactions and 
probabilities depending on arrangement of three ions in 
Table ini From Tabic El the probabilities of FM and AF 
effective interactions are 1 — 2a;(l — x)^ and 2x(l — x)^, 
respectively. The spatial average of signs of effective in- 
teractions between NNLs, cr""\ is given by 

^ = 4a;^ _ Sx^ +4x - 1. (5) 

From this equation, it is clear that ct""' is negative re- 
gardless ofx(0<x<l). Thus, the effective interaction 
between NNLs is always FM. From this fact, it follows 
that the spin arrangement between odd-numbered (even- 
numbered) layers is FM along the interlayer direction in 
ordered phases. Note that the ferromagnetic correlation 
between NNLs exists in not only collinear spin structures 
but also non-collinear spin structures. This is because 
thermal and disorder fluctuations cause the ferromag- 
netic effective interaction regardless of spin structure. 

To compare the properties of site-random model and 
bond-random (BR) model where FM and AF interactions 



TABLE L Decision rule of interactions for combinations of 
ions (J > 0). 





intralayer 


interlayer 


A-A 


+J(AF) 


+J(AF) 


A-B 


+JiAF) 


-J(FM) 


B-B 


+J(AF) 


-J(FM) 



TABLE II: Effective interactions and probabilities of next 
nearest-neighbor spin pairs along the interlayer direction de- 
pending on the arrangement of three ions. 



arrangement 


effective interaction 


probability 


A-A-A 


FM (a?,"' 


= -1) 


{l-xf 


A-A-B 


AF {afi' 


= +1) 


x{l - xf 


A-B-A 


FM (a-r' 


= -1) 


x{l ~ xf 


B-A-A 


AF (a-\"' 


= +1) 


x{l - xf 


A-B-B 


FM {af^' 


= -1) 


x^{l - x) 


B-A-B 


FM 


= -1) 


x^{l - x) 


B-B-A 


FM (a'^ 


= -1) 


x^{l - x) 


B-B-B 


FM (a-^' 


= -1) 


x^ 



are placed randomly on each bond, we analyze the effec- 
tive interactions for bond-random model in a similar way. 
The spatial average of signs of interlayer interactions CTbr 
and that of effective interaction between NNLs cr™' are 
given by 

0^=1-2?/, (6) 
^' = Ay^-Ay + l, (7) 

where y is the concentration of FM bonds along the in- 
terlayer direction, and is at y = y* := 1/2. Since 

a^^^ is equal to a.t y — y* , the bond-random model does 
not have long-range order in the interlayer direction at 
y — y* . Thus, the existence of the FM effective inter- 
action between NNLs regardless of x is a characteristic 
feature of the site-random model and is crucial for ob- 
taining a finite transition temperature even at x = x* in 
our model. 



IV. SIMULATION RESULTS 

We study the finite temperature properties of the clas- 
sical Heisenberg model with site-random interlayer cou- 
plings defined by Eq. (IT])onaniV = _LxLxL simple 
cubic lattice with a periodic boundary condition. We use 
Monte Carlo simulations in which the standard heat-bath 
method is adopted. Before starting the simulations, we 
specify the B ion concentration x such that the number 
of B ions (A^B iV x x) is an integer. In each simulation, 
we prepare random configurations of A and B ions and 
set the interaction of each bond according to Table H] 
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FIG. 5: (Color online) Phase diagram of temperature T/J 
versus B ion concentration x obtained from Monte Carlo sim- 
ulations. Asterisk and arrow denote parameters used in Sec. 
IV B (see Fig. |6]) and C (see Fig. llOp . respectively. In the 
mixed phase, (nnn) and (ttttO) orders coexist, and the spin- 
ordering pattern is the "random fan-out state" as shown in 
Fig.El 



FIG. 6: Correlation function between two spins along the 
interlayer direction (c-axis) for L = 24 at T/J = 0.1 obtained 
from Monte Carlo simulations. Dotted line is visual guide. 



B. Spin-ordering pattern in mixed phase 



A. Phase diagram 

In this section, we draw the phase diagram of tem- 
perature T/J versus B ion concentration x. To deter- 
mine transition temperatures, we calculate the Binder 
ratio U4{q) of magnetization vector m{q) characterized 
by wave vector q: 



1 



fniq) — > s^e" 



N 



UA{q) 



m{q)\' 



\m{q) 



(8) 



(9) 



Here, (• • ■ ) is the thermal average and [■ ■ • ]av is the ran- 
dom average over arrangements of ions. Hereafter, 
represents the position vector of the site i. We prepare 
64-1024 random configurations of ions. The transition 
temperature is determined at the crossing point of Ui{q) 
by using L = 12-32 data. The phase diagram obtained 
from Monte Carlo simulations is shown in Fig. [5l There 
are three ordered phases in the phase diagram. Here, we 
determine each phase structure according to the struc- 
ture factor: 



S{q):^N {\m{q)\ 



(10) 



In the (tttttt) and (vrTrO) ordered phases, the structure 
factor has peaks only at each corresponding wave vector. 
In the mixed phase, structure factor peaks at two distinct 
wave vectors (tttttt) and (ttttO) develop. There are two 
phase boundaries in the phase diagram. The blue circles 
and red pentagons in Fig. [S] are determined by Ui{mnr) 
and Ui{7nrO), respectively. 



We next investigate the spin-ordering pattern in the 
mixed phase. In the mixed phase, peaks of structure 
factors at (tttttt) and (ttttO) develop; namely, Neel or- 
der appears in each layer. To clarify the spin structure 
along the interlayer direction, we calculate a correlation 
function between two spins along the interlayer direction 
(c-axis), which is defined by 



1 

N 



^ {s{ri) ■ s{ri + rcBc)) 



(11) 



Here, Cc is a unit vector in the interlayer direction and 
Tc is the distance between two spins. Furthermore, s{ri) 
denotes the spin at position r^, which is the same as s^. 
We adopt x = 19/64 = 0.296875, which is slightly larger 
than X* ~ 0.2929. Figure [5] shows the dependence of 
distance of G"'(rc) for L = 24 at T/J = 0.1 which is 
below the transition temperature. The spin correlation 
along the interlayer direction is FM for an even number of 
Tc. This result is consistent with the expectation based 
on the effective interaction between NNLs discussed in 
Sec. HI. In contrast, the spin correlation is nearly zero 
for an odd number of rz- This means that the angle be- 
tween nearest-neighbor spin pairs along the interlayer di- 
rection is nearly tt/2. The value of spin correlation for an 
odd number of r,. corresponds to the angle between spins 
belonging to odd-numbered and even-numbered layers. 
Next, we consider the x dependence of this angle. 

We calculate the x dependence of the angle 9 between 
the staggered magnetization vector of an odd-numbered 
layer rrio and that of an even-numbered layer rUe in the 
ground state. The magnetization vectors mo, mo, and 
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FIG. 7: (Color online) (a) Dependence on concentration 
X of the angle 9/n in the ground state, (b) Dependence 
on concentration x of magnitude of magnetization rriav := 
(|mo| + |mc|)/2 in the ground state. 
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where q = (vrvrO). We obtain the spin configuration in 
the ground state by the steepest descent method corre- 
sponding to zero temperature simulation. We prepare 
64 random configurations of ions and use the snapshot 
obtained from Monte Carlo simulations aXT/ J — 0.1 as 
the initial spin configuration. Figure [7] (a) shows the x 
dependence of the angle 0/tt \x\ the ground state. As 
the value of x increases, the value of 9 /i: changes from 1 
corresponding to the (tttttt) order to corresponding to 
the (TTTrO) order. Figure [7] (b) shows the x dependence of 
TTiav := (|"T'o| + I'TT'cD/Z in the ground state. The value 
of rriav is less than 1 in the mixed phase. This means 
that the spin arrangement of each layer is not perfect 
Neel order. Spin directions are distributed in a fan-shape 
around each magnetization axis owing to the effect of 
random fields from neighboring layers. This fan-shaped 
distribution of spins can be seen in the snapshot of spin 
directions (Fig. [5]), which are projections onto the plane 
perpendicular to the vector mo x rric- Thus, jmol and 
I mo I indicate the shape of the spreading spin- fan. Note 
that the absolute values of mo and rric should be the 
same, |mo| ~ |mo|, because of the symmetry. This fact 



FIG. 8: (Color online) Snapshots of spin directions for x = 
0.296875 and 0.53125 when the lattice size is L = 12 in the 
ground state. Coordinates (sx, Sy) define the orthogonal plane 
of the vector rUo x me. 
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Random Fan-Out State 



FIG. 9: (Color online) (a) "Average" spin directions in the 
spin configuration of the random fan-out state. In each layer 
(afo-plane), Neel order appears. Along the interlayer direction 
(c-axis), the angle 9 between nearest-neighbor "average" spin 
pairs changes from tt to with increasing x as shown in Fig. [7| 
(a). The correlation between NNLs is FM. (b) Individual 
spins are randomly directed around the average direction. 



is adopted in the mean-field analysis as an assumption, 
which is discussed in Sec. V. 

From these results, the schematic spin configuration in 
the mixed phase can be depicted as shown in Fig. [HI Fig- 
ure [2] (a) shows a schematic configuration of spin axes mo 
and mo. The direction of each arrow is the actual "aver- 
age" magnetic direction. In this configuration, the spins 
are grouped into four sublattices according to whether 
they belong to an even/odd layer and which of the two 
sublattices in each layer. Each sublattice corresponds to 
a distinct average spin direction as shown in Fig. |9] (a) . 
Individual spins spread in the shape of fan around the 
average direction as shown in Fig.[H](b). From these fea- 
tures, we refer to this bulk spin-ordering pattern in the 
mixed phase as the "random fan-out state" . This random 
fan-out state is an exotic bulk spin ordering that explains 
the simultaneous appearance of (tttttt) and (ttttO) wave 
vectors without any phase separation. In this model. 
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(b) lower-temperature transition 
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FIG. 10: (Color online) (a) Finite-size scaling plots of the higher-temperature phase transition from the paramagnetic phase 
to the (tttttt) ordered phase. We adopt the critical exponents of three-dimensional Heisenberg universality class {i/ = 0.704, 
rj — 0.025). (b) Finite-size scaling plots of the lower-temperature phase transition from the (tttttt) ordered phase to the mixed 
phase. We adopt the critical exponents of three-dimensional XY universality class {f = 0.672, rj = 0.038). 



the ferromagnetic correlation between NNLs exists as ex- 
plained in Sec. III. 



C. Universality classes of phase transitions 

We study the universality classes of phase transitions 
of our model. In the phase diagram (Fig.[5|), there are two 
types of phase boundaries. To make clear the universality 
classes of each phase transition, x is set to 3/16 = 0.1875 
such that transition temperatures are separated suffi- 
ciently. For this parameter, the intermediate phase is 
the (tttttt) ordered phase (see the dotted arrow in Fig. [5]). 

First, we investigate the higher-temperature phase 
transition from the paramagnetic phase to the (tttttt) or- 
dered phase. From the Harris criterion (36j. we expect 
that the higher-temperature phase transition belongs to 
the three-dimensional Heisenberg universality class. This 
is because the critical exponent a is negative in the 
three-dimensional Heisenberg model, and thus the dis- 
order should not affect the universality class. To obtain 
the transition temperature and confirm the critical expo- 
nents, we calculate the correlation function G'^(rc), which 
is defined by Eq. pT|) . The finite-size scaling relations of 
(i-dimensional systems are given by 



G%L/2) 



G"=(i/4) 



(15) 
(16) 



where ^ and $ are scaling functions and t := T — Tf |37|. 
We determine the transition temperature Tc as the cross- 
ing point of G"=(L/2)/G"=(L/4) using L = 20-32 data 
and obtain T^/J = 1.085(5). This transition temper- 
ature is consistent with the one found for J74 (tttttt) in 



Sec. IV A. The finite-size scaling using the critical ex- 
ponents of the three-dimensional Heisenberg universality 
class {v = 0.704, rj = 0.025) J3| are shown in Fig. [TO] (a). 
Since the data are well fitted by scaling relations, we 
conclude that the higher-temperature phase transition 
belongs to the three-dimensional Heisenberg universality 
class in accord with the Harris criterion. 

Next, we investigate the lower-temperature phase tran- 
sition from the (tttttt) ordered phase to the mixed phase. 
The (tttttt) ordered phase is translationally symmetric 
with the 0(3) spin rotation symmetry broken down the 
spin rotation symmetry U(l). In the mixed phase, both 
the translational symmetry and the U(l) spin rotation 
symmetry are broken. Thus, we expect that the transi- 
tion to the mixed phase is characterized by the breaking 
of U(l) spin rotation symmetry. In other words, we ex- 
pect that the lower-temperature phase transition belongs 
to the three-dimensional XY universality class. To obtain 
the transition temperature and confirm the critical ex- 
ponents, we calculate the magnetization vector m(TTTTO) 
defined by Eq. ([8]). The finite-size scaling relations are 
given by 



]m(q)P)2 



;,(g)^A.i^!iMQoci-".g(tLV.), 



T 



(17) 



(18) 



where / and g are scaling functions, and q = (ttttO). 
We determine the transition temperature as the cross- 
ing point of C/4(ttttO) using L = 20-32 data and obtain 
Tc/J = 0.4585(5). The finite-size scaling using the crit- 
ical exponents of the three-dimensional XY universality 
class (i/ = 0.672, ?] = 0.038) @ are shown in Fig. [TO] 
(b). Although we obtain a reasonably good fit, it is 
not good enough to detect the small difference between 
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XY critical exponents and the other critical exponents 
in a three-dimensional system. However, from the view- 
point of spin rotation symmetry, we can deduce that the 
lower-temperature phase transition belongs to the three- 
dimensional XY universality class. 



V. MEAN-FIELD CALCULATIONS 



component h in the definition of hi, because h does 
not contribute to the free energy. 

The excess free energy of a one-layer system from the 
effects of molecular fields hi is given by 

'^n^) = -^EE*^^x^^-'^^-E^" -"^^ (24) 

i j i 



In this section, to obtain an intuitive understanding of 
the emergence mechanism of the mixed phase, we inves- 
tigate the effect of random interlayer couplings by mean- 
field calculations. For simplicity of notation, we study 
the system where the intralayer interactions are FM. Un- 
der the gauge transformation at alternating sites, this 
model is equivalent to our model given by Eq. ([1]) in the 
case of no external field. By applying the inverse of the 
gauge transformation to this model, we can obtain the 
same results as from the original model. In this section, 
we define rrio as the uniform magnetization vector of an 
odd-numbered layer and rrio as that of an even-numbered 
layer. The magnitudes of these magnetization vectors are 
assumed to be the same (m := \mo\ = Imd), which is a 
reasonable assumption as stated in Sec. IV B. Here, let 
9 be the angle between mo and rrie. In this mean- field 
calculation, to consider a simpler explanation, we focus 
on one odd-numbered layer with molecular fields from 
neighboring layers. This is because the circumstances 
of odd-numbered layer and even-numbered layer are the 
same, and thus it is enough to investigate either layer. 
The molecular field from neighboring layers at the site i 
in an odd-numbered layer is given by 



h, 



-2 J'tJiTne 



(19) 



where J' (> 0) denotes the nearest-neighbor interaction 
along the interlayer direction, and cr^ denotes the sign of 
interaction with neighboring interlayer sites. Note that 
we have replaced the spins on the neighboring layers with 
the average spins me and neglected the fluctuation, and 
the contributions of upper and lower layers are assumed 
to be the same. Here, we separate ai into the spatial 
average a and the deviation Adi from a depending on 
the site {ai — a + Aai). The random average of deviation 
[ACTiJav = is assumed. Accordingly, hi can be separated 
into 



h, = ft," + Aftf 



Ahf 



where 



h" 
Aftj' 
Ahi 

,11 



-2J'mcos9ael, 
-2J'm cos 0Aaiel , 
-2J'm sin 9A<7ie;^. 



(20) 

(21) 
(22) 
(23) 



The symbols ei and e^, respectively, are unit vectors 
that are parallel and perpendicular direction of mo to 
mo. Here, we did not include the uniform transverse 



where i and j denote site indexes on one layer, and the 
susceptibility tensor is defined as 



drrii 



(25) 



and nrii is the magnetization vector at the site i. The 
random average of the product of deviations is expressed 
as 



[Aa.AajU = (1 - ^')% =: A(t2%. (26) 

Hence, we obtain the free energy per site by taking the 
random average as following expression: 



Sf := [SF{0)U/N 



2 3 



Ah\ 



i 

T'\ II -2 

1 -^uniform ^ 



^ ft" • rrii 



N 



(xLal - Xlical) A(72|cOs2 



(27) 



where N is the number of sites in one layer. The symbols 

xllniform' xical', and xi^cai represent the uniform longi- 
tudinal susceptibility, the local longitudinal susceptibil- 
ity, and the local transverse susceptibility, respectively: 

Xuniform SjXij; Xlocal Xii i ^ud Xlocal Xii 

all i. In the mean-field approximation, the magnitude of 
magnetization m in the field h along the magnetization 
vector is given by the Langevin function: 



L(ph) = coth(/3/i) - ^ ^m. 

ph 



(28) 



The molecular field from the same layer is given by 

h^^^ = 4Jm. (29) 
Thus, the susceptibilities obtained by mean-field approx- 
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imation are calculated by using 



^uniform 



-^iocal 



1 - {h^^lm)^l 



cal 



II _ dLjph) 

Xlocal 



h=h^ 



Xlocal 



(30) 

(31) 
(32) 



We calculate the angle that minimizes the excess free 
energy expressed by Eq. (P?]) . The magnitude of the mag- 
netization TO is obtained from the self-consistent equa- 
tion: 



L {Ph 



(33) 



The phase diagram of temperature T / J versus B ion con- 
centration X for J = J' is shown in Fig. [TT] (a). In this 
phase diagram, there is a phase corresponding to the 
mixed phase where tMo and rrif. are not parallel or an- 
tiparallel to each other. We use the molecular field from 
the same layer only when we calculate the magnitude of 
magnetization m. Then, in the mean-field calculations, 
the transition temperatures of the higher-temperature 
transition are clearly the same regardless of x. To ex- 
amine the properties of this phase, we calculate the x 
dependence of 9/n in the ground state (Fig. [TT] (b)). As 
the value of x increases, the value of 9/-K changes from 
1 corresponding to the (tttttt) order to corresponding 
to the (ttttO) order. Since the phase diagram and the 
behavior of 9 obtained from mean-field calculations are 
qualitatively consistent with the Monte Carlo simulation 
results, we expect that the emergence mechanism of the 
mixed phase is explained by mean-field calculations. 
Next, we consider the stabilization mechanism of the 

mixed phase. From Eq. ([27]), when X\oca\ > xlLcai ^'^'^ 
(Aft,^)^ > 0, the excess free energy is minimized when 
6'/7r 7^ 1 or 0. Thus, we conclude that the mixed phase 
emerges by the following mechanism, (i) In the pres- 
ence of finite magnetization, spins are more susceptible 
to transverse field than longitudinal field, which is physi- 
cally natural, (ii) Transverse fields from neighboring lay- 



ers exist as random fields Aft, 



which come from the fluc- 
This condition is 



tuation of random interactions Act^ 
realized by the existence of random interlayer couplings. 
Thus, we clarify that the mixed ordering, in which the 
magnetization vectors of neighboring layers are not par- 
allel or antiparallel, is induced by random interlayer cou- 
pling. 



VI. DISCUSSION 




FIG. 11: (Color online) (a) Phase diagram of temperature 
T / J versus B ion concentration x for J = J' obtained from 
mean-field calculations, (b) Dependence on concentration x 
of angle 6/ti between ruo and rric in the ground state. The 
points denote the simulation results for L = 12 and 24, which 
are the same results shown in Fig. [7] (a). 



configuration obtained from Monte Carlo simulations at 
low temperature. 



A. Robustness of random fan-out state 

In this subsection, we discuss the relationship between 
decision rules of interactions and the spin-ordering pat- 
tern. In general, the spatial average of interlayer random 
interactions J and its fluctuation A J corresponding to 
the random fields depend on the signs and absolute val- 
ues of interactions: 



J ~ ^aa(1 — ■ 



(34) 



A J = \/ JiA(l - + 2Jl^x{\ -x) + Jl^x'^ - J 



-2 

(35) 



where Jaa, Jab, and Jbb denote the interactions be- 
tween A- A, A-B, and B-B along the interlayer direction, 
respectively. Then, the spatial average of effective inter- 
action between the NNLs, J""', is given by 



J™i cx - JXa(1 - xf - 2JaaJabx{1 - x)' 
- Jabx{^ - ^) ~- 2Jab>/bb2^^(1 - x) 



'BB-^ 



(36) 



In this section, we examine whether the random fan- 
out state is robust with respect to the decision rules of 
interactions and anisotropy. We also show the results 
of Rietveld analysis of Sr(Feo.7Mno.3)02 using the spin 



Here, we assume that Jaa is always AF and the absolute 
values of interactions are the same or zero. Then, four 
rules in Table Hill are all possible decision rules of interac- 
tions where FM and AF interactions exist and J""' < 
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TABLE III: Decision rules of interactions along the interlayer direction. Rule 1 is adopted in Sec. Ill, IV, and V. 





Jaa 


Jab 


Jbb 


J/J 


AJ/J 


J^/J oc 


X* 


Rule 1 


+J(AF) 


-J(FM) 


-J(FM) 


2a;^ ~4x + l 




4a;(l -xf -I 


(2 - x/2)/2 


Rule 2 


+J(AF) 


-J(FM) 





- 4x + 1 


{l-x^-T)i 


- + 4a; - 1 


1/3 


Rule 3 


+J(AF) 


+J(AF) 


-J(FM) 


-2x-^ + 1 




-4x^ + 4x^-1 


l/^/2 


Rule 4 


+J(AF) 





-J(FM) 


-2x + 1 


{l-2x + 2x' - f)^ 


-3a;^ + 3a; - 1 


1/2 



regardless of x. Note that rule 1 is adopted in the pre- 
vious sections. Under these rules and using Monte Carlo 
simulations and mean-field calculations, we find the ap- 
pearance of the mixed phase, where the spin-ordering 
pattern is the random fan-out state. Furthermore, we 
also find that (tttttt) and (ttttO) magnetic peaks develop 
at the same temperature and the angle between the stag- 
gered magnetization vectors of neighboring layers is 7r/2 
at a: = X* for all rules shown in Table IIIII 

Next, we consider the effect of easy-plane anisotropy. 
We study the finite temperature properties of the XY 
spin version of our model by using Monte Carlo simu- 
lations. Using rule 1 in Table HID we find the mixed 
phase where spin-ordering pattern is the random fan-out 
state. Thus, we conclude that the random fan-out state 
is not prohibited by the effects of easy-plane anisotropy. 
In contrast, the random fan-out state is prohibited by 
the effects of easy-axis anisotropy, because the energy is 
raised when rrio and mo are not parallel or antiparallel 
to each other. 

From these discussions, the random fan-out state is 
robust with respect to the details of our model, which 
satisfies the following conditions, (i) The random FM-AF 
interactions exist along the interlayer direction, (ii) The 
effective FM interaction between NNLs exists regardless 
of X. (iii) There is no easy-axis anisotropy. 

B. Relation with experiments 

We discuss the spin structure in the low-temperature 
state of Sr(Feo.7Mno.3)02. In neutron diffraction mea- 
surements, magnetic peaks at (tttttt) and (ttttO) develop 
at low temperatures as stated in Sec. II. This result can 
be explained by two possible scenarios — phase separa- 
tion or a bulk spin structure that is the random fan-out 
state. In the former case, from the ratio of intensities 
at (tttttt) and (ttttO), we can determine the volume frac- 
tion of each phase. In the latter case, the angle 9 be- 
tween the staggered magnetization vectors of neighbor- 
ing layers can be estimated from the ratio of intensities 
at (tttttt) and (ttttO). Therefore, an important contrast 
between the two scenarios lies in the ordering tempera- 
ture. In the phase separation scenario, ordering of the 
two magnetic phases should occur at different temper- 
atures because the absolute values of each interaction 
of Sr(Feo.7Mno.3)02 should be different. On the other 
hand, in the bulk spin scenario, even when the absolute 



values of each interaction are different, there should be 
a cross point of the phase boundary between the (tttttt) 
and (ttttO) orders in the phase diagram. Based on ex- 
perimental data. Fig. |4] shows that the magnetic peaks 
at (tttttt) and (ttttO) lose their intensities at nearly the 
same temperature, which is a strong indication that the 
low-temperature state of Sr(Feo.7Mno.3)02 is likely to be 
the random fan-out state. At the cross point, the angle 
between the staggered magnetization vectors of neigh- 
boring layers should be nearly tt/2 and the two magnetic 
peaks should disappear at nearly the same temperature. 
Thus, if the angle 9 obtained by Rietveld structural re- 
finement is nearly tt/2, the simultaneous emergence of the 
two types of wave vectors can be naturally explained by 
the bulk spin structure scenario. To determine the angle 
in the random fan-out state, we perform Rietveld struc- 
tural refinement of Sr(Feo.7Mno.3)02 using the random 
fan-out state. The results of Rietveld structural refine- 
ment of the neutron diffraction data and full details of 
the refined parameters are shown in Fig.[T2]and Table HVl 
respectively. The neutron diffraction pattern is well fit- 
ted by using the random fan-out state with 9 = 84° in 
Fig. [121 Thus, we conclude that it is highly possible that 
the mixed ordering in Sr(Feo.7Mno.3)02 is due to the bulk 
spin structure, which we call the random fan-out state. 

Since the spin-ordering pattern in Sr(Feo.7Mno.3)02 at 
low temperature is explained by our model, we expect 
that the x — T phase diagram of Sr(Fei_xMna;)02 will 
agree with the theoretical phase diagram (Fig. [5]) qualita- 
tively. Indeed, it appears that the transition temperature 
from the paramagnetic phase to the (tttttt) ordered phase 
Tn decreases with increasing x when x = 0.0, 0.1, and 
0.2 as stated in Sec. II. This behavior of Tn is consistent 
with the theoretical phase diagram (Fig. [5])- If this theo- 
retical phase diagram is valid for Sr(Fei_2:Mnj:)02, it will 
be observed that the (ttttO) peaks in Sr(Fei_xMn2:)02 for 
X = 0.1 and 0.2 develop at low temperatures. In future 
work, we will perform neutron diffraction measurements 
on Sr(Fei_a;Mna;)02 for several x in addition to cc = 0.3 
at low temperature. Moreover, we will study the rela- 
tionship with the theoretical interpretation discussed in 
this paper. 

VII. SUMMARY AND CONCLUSION 

This paper is summarized as follows. In Sec. II, 
we have reported the neutron diffraction pattern of 
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FIG. 12: (Color online) Observed (crosses), calculated (lines), 
and difference plots from the structural refinements of powder 
neutron diffraction data for Sr(Feo.7Mno.3)02 at 10 K. The 
calculated data (lines) denote the intensity of the random fan- 
out state. The upper and lower ticks represent the positions of 
the calculated chemical and magnetic reflections, respectively. 



TABLE IV: Rietveld refinement for Sr(Feo.7Mno.3)Oi. 



10 K 



RT" 



a (A) 
c(A) 

B in Sr (A) 

B in (Feo.rMno.s) 

B in O (A) 

Moment {/j-b) 

Angle 9 

-Rwp (%) 

Rp (%) 

,,2 



(A) 



4.0014(1) 

3.4406(1) 

1.31(4) 

1.08(4) 

1.51(5) 

2.0 

84° 

15.6 % 
15.2 % 
3.85 



4.0055(1) 

3.4553(1) 

0.42(7) 

0.21(6) 

0.91(5) 





9.11 % 
7.01 % 
3.77 



" All the refinements were performed by using the PA/mmra 
space group with Sr on Id (1/2,1/2,1/2), (Feo.rMno.s) on la 
(0, 0, 0), and O on 2f (1/2, 0, 0). 
The refinement for RT is from Ref. 



Sr(Feo.7Mno.3)02. We have observed that magnetic 
peaks corresponding to (tttttt) and (ttttO) wave vectors 
emerge simultaneously at Tn — 240 K. 

In Sec. Ill, to investigate whether there is a bulk spin 
structure that explains the simultaneous appearance of 
(tttttt) and (ttttO) wave vectors, we have introduced the 
classical Heisenberg model with site-random interlayer 
couplings on the cubic lattice as a simplified effective 
model of Sr(Fei_a;Mna;)02. This model consists of two 
types of magnetic ions, labeled A and B, that correspond 
to Fe and Mn ions. The interactions in a plane have been 
assumed to be uniform antiferromagnetic. The interac- 
tions along the interlayer direction between A ions have 
been assumed to be antiferromagnetic, and other inter- 



actions have been assumed to be ferromagnetic. In this 
decision rule of interactions, the ferromagnetic effective 
interaction between next-nearest layers along the inter- 
layer direction exists regardless of the B ion concentra- 
tion X. This is a characteristic feature of site-random 
interlayer couplings. 

In Sec IV, by performing Monte Carlo simulations, we 
have shown the existence of a bulk spin structure that ex- 
plains the simultaneous appearance of (tttttt) and (ttttO) 
wave vectors without any phase separation. In this spin 
structure called the random fan-out state, there are three 
features: (i) the arrangement of spins in each layers is 
Neel order, (ii) correlation between next-nearest layers 
has ferromagnetic, (iii) the magnetization vector in odd- 
numbered layers and that in even-numbered layers are 
not parallel or antiparallel to each other. As the value 
of X increases, the angle between these ordering vectors 
changes from tt corresponding to the (tttttt) order to 
corresponding to the (ttttO) order. Moreover, we have 
constructed the phase diagram of temperature versus B 
ion concentration. In the phase diagram, the succes- 
sive phase transitions have been observed. We have also 
found the critical exponents of each phase transition. 

In Sec. V, to clarify the emergence mechanism of the 
mixed phase, we have investigated the effects of random 
interlayer couplings by using mean- field calculations. It is 
important for the existence of a mixed phase that spins 
be more susceptible to transverse field than longitudi- 
nal field. Since the transverse component of the random 
field from neighboring layers becomes large in the mixed 
phase, the magnetization vectors of neighboring layers 
are not parallel or antiparallel. 

In Sec. VI, we have discussed whether the random 
fan-out state is robust with respect to the details of 
the model. We have clarified that the random fan- 
out state becomes stable, when random ferromagnetic- 
antiferromagnetic interactions exist and the ferromag- 
netic effective interaction between next-nearest layers ex- 
ists regardless of x. Furthermore, we have determined 
that the random fan-out state is not prohibited by the 
effects of easy-plane anisotropy. We have also performed 
Rietveld structural refinement of Sr(Feo.7Mno.3)02. The 
neutron diffraction pattern is well fitted by using the 
random fan-out state. Thus, we conclude that the ran- 
dom fan-out state is a reasonable spin-ordering pattern 
of Sr(Feo.7Mno.3)02 at low temperature. 

In this research, we have found a bulk spin structure 
called the random fan-out state that explains the simulta- 
neous appearance of (tttttt) and (ttttO) wave vectors with- 
out any phase separation. The model has been intro- 
duced as a simplified effective model of Sr(Fei_3;Mn3;)02, 
and the neutron diffraction pattern of Sr(Feo.7Mno.3)02 
is well fitted by using this spin structure. In our model, 
we have introduced the decision rule of interactions which 
induces disorder of interaction along only the interlayer 
direction. Thus, the spatial distribution of frustration in 
this model is different from that in systems having an 
isotropic crystal structure. Indeed, in a plane, since frus- 
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tration does not exist, Neel order appears. In contrast, 
random ferromagnetic-antiferromagnetic interactions ex- 
ist along the interlayer direction. Thus, effect of frustra- 
tion appears only between layers. At low temperature, 
the spin structure appears in which the angle between 
the neighboring layer magnetization axes are the same. 
This structure is reminiscent of the spiral spin structure 
induced by the effects of frustration. However, in our 
site-random model, the interactions along the interlayer 
direction are correlated, and ferromagnetic correlation 
between next-nearest layers appears. Thus, a spiral spin 
structure characterized by a single wave vector does not 
appear, but an unusual spin-ordering pattern expressed 
by two wave vectors appears. Furthermore, in our model, 
the spin-fan structure exists in the ground state without 
an external field. In frustrated magnetic systems, the 
well-known spin-fan structure is realized by applying an 
external field along the direction parallel with the spiral 
plane 40 - 42] . However, since the spin- fan structure of our 
model arises from an intrinsic effect of site-random cou- 
plings, the emergence mechanism of the spin- fan struc- 
ture is different from that of the conventional one (40l - l42j . 
In our model, by the fluctuation of random interactions, 
the random fields from neighboring layers exist depend- 
ing on the site, and thus spin directions have a fan-shaped 
distribution around each magnetization axis. 

The random fan-out state is a novel type of spin- 
ordering pattern that comes from site-random interlayer 
couplings. Thus, it is important to investigate the re- 
sponse to magnetic field or other external fields and spin 
wave excitation of this spin structure. Recently, many ex- 
otic physical properties induced by competition and har- 
monization between spin structure and other degrees of 
freedom have been report ed jisl - iisj . and thus the search 
for new physical properties hidden in the random fan- 



out state is extremely important from the perspective of 
materials science. 

In statistical physics, the phase transition and char- 
acteristic spin structure in anisotropic random systems 
will also be important areas of investigation. Although 
dynamical properties in random-spin systems have been 
studied for a long time, the relation between the dynam- 
ics and spatial distribution of disorder has not been estab- 
lished yet to the best of our knowledge, and thus remains 
a challenging problem to be investigated. 
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